2D skew scattering in the vicinity and away from resonant scattering condition 
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We studied the energy dependence of the 2D skew scattering from strong potential, for which 
the Born approximation is not applicable. Since the skew scattering cross section is zero both 
at low and at high energies, it exhibits a maximum as a function of energy of incident electron. 
We found analytically the shape of the maximum for an exactly solvable model of circular-barrier 
potential. Within a rescaling factor, this shape is universal for strong potentials. If the repulsive 
potential has an attractive core, the discrete levels of the core become quasilocal due to degeneracy 
with continuum. For energy of incident electron close to the quasilocal state with zero angular 
momentum, the enhancement of the net cross section is accompanied by resonant enhancement of 
the skew scattering. By contrast, near the resonance with quasilocal states having momenta ±1, the 
skew scattering cross section is an odd function of energy deviation from the resonance, and passes 
through zero, i.e., it exhibits a sign reversal. In the latter case, in the presence of the Fermi sea, the 
Kondo resonance manifests itself in strong temperature dependence of the skew scattering. 

PACS numbers: 72.25.Dc, 73.43.-f, 72.25.Rb, 72.20.Dp 
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I. INTRODUCTION 

The phenomenon of skew scattering has its origin in 
a spin-orbit part oc er • [k x VV(p)] of the scattering po- 
tential, V(p). Due to this part, the scattering ampli- 
tude between the states with momenta k and k' acquires 
a contribution, f a oc <x • [k x k'], which is asymmetric 
with respect to the scattering angle, 9, between k and 
k'. Therefore, in two dimensions, differential cross sec- 
tion acquires an asymmetric contribution 



da a 
~~aW 



25fte[(/ s )*/ a ] oc asin( 



(1) 



where f s is the symmetric part of the scattering ampli- 
tude, i.e., the scattering amplitude in absence of spin- 
orbit potential, and a = ±1 is the spin projection on the 
normal to the 2D plane. 

Skew scattering, introduced almost 80 years ago^, 
had recently attracted a lot of interest, since it is a key 
ingredient of the anomalous Hall e&ect£££&£& as well 
as of the spin-Hall eSec^dL^W^W 7 , Detailed 
theoretical calculations of a a were carried out for atomic 
systems! 8 . Complexity of these calculations stems from 
the fact that a a = in the lowest Born approximation 
for f s . The results of the second Born approximation for 
f s , yielding a finite a a , are presented in Ref. |l8| for the 
model of the Thomas-Fermi screening of the charge of a 
nucleus. 

Calculations of Ref. [Hwere recently utilized in Ref. 
to estimate the magnitude of skew scattering by a donor 
impurity in 2D electron gas. Concerning the energy de- 
pendence of the skew-scattering, a a (E), in calculation of 
the anomalous Hall and spin-Hall effects, it should be 
taken at E — Ep, where Ep is the Fermi energy. Thus, 
the skew scattering is sensitive to the electron density 
via Ep. Then the question about the explicit form of 
a a (E) in two dimensions arises. Nontriviality of this de- 
pendence stems from the fact that at low energy f s (E), 
is dominated by the angular momentum, 1 = 0. As a 
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FIG. 1: Schematic plot of circular-barrier potential with 
height, Vo, and radius, b. Attractive core at p < a contains a 
localized state with energy, Eo, which has a width, V, due to 
degeneracy with continuum. 



result, the skew-scattering cross section, a a (E), turns to 
zeroi at E — > 0. Since a a (E) turns to zero at large E as 
well, we conclude that it should pass through a maximum 
at a certain finite energy. 

Study of the energy dependence of the skew scattering 
is the focus of the present paper. We demonstrate that 
this dependence exhibits especially rich behavior for a 
"strong" scattering potential with characteristic magni- 
tude, Vo, and characteristic radius, b, satisfying the con- 
dition Vo 3> h 2 /mb 2 , where m is the electron mass. For 
such potentials, the Born approximation does not apply 
at low energies. Then the calculation of the skew scat- 
tering requires the knowledge of scattering phases in the 
absence of spin-orbit coupling, while the spin-orbit cou- 
pling should still be treated perturbatively. 

Even more interesting energy dependence of the skew 
scattering emerges in the case when the scattering po- 



2 



tential possesses a core, as shown in Fig. 1. Then the 
quasilocal state, E , due to this core, has a small width, 
r. It is known, that spin-independent cross section ex- 
hibits a resonant behavior in the vicinity of Eq. We will 
demonstrate that skew scattering is also resonantly en- 
hanced near E — Eq. A particular interesting situation 
realizes when the quasilocal state at E = E\ corresponds 
to nonzero momenta, I = ±1. The peculiarity of this 
situation is that weak spin-orbit-induced splitting of the 
states with / = 1 and / = — 1 can be comparable to T. 
Then the dependence, a a (E), is nonperurbative with re- 
spect to spin-orbit coupling. Moreover, we demonstrate 
that the skew scattering can change sign within a narrow 
energy interval around E = E\. Finally we show, in the 
presence of the Fermi sea, the Kondo physics, that mod- 
ifies the elastic scattering with / = also modifies the 
energy dependence of the skew scattering near the Fermi 
level within the narrow interval (E — Ep) of the order of 
the Kondo temperature. 

To the best of our knowledge, the "exact" (without 
utilizing the Born approximation) energy dependence of 
the 2D skew scattering has been previously studied by 
Hankiewicz and Vignalei^ and by Huang et al~ In both 
papers the study of a a (E) was carried out without ex- 
pansion with respect to a weak spin-orbit coupling; the 
results for the skew scattering are presented in the form 
of numerical plots for certain parameters of the scatter- 
ing potential and the magnitude of the spin-orbit cou- 
pling. This complicates the inference of general features 
of a a (E). We, on the other hand, made use of the small- 
ness of the spin-orbit coupling at the first step. This 
allowed us to obtain closed analytical results for the skew 
scattering as a function of energy. For strong scatter- 
ing potential, V(p), we show that the energy dependence 
of <r a is universal within two parametrically wide do- 
mains of low energies before the Born approximation be- 
comes valid. Characteristics of V(p) determine the en- 
ergy scales, but not the form of a a (E) in these domains. 
Concerning the resonances, around energies E — Eq and 
E = Ei of the quasilocal states, a a (E) exhibits a uni- 
versal behavior as well, in the same way as the spin- 
independent scattering. Since quasilocal states manifest 
themselves only in the narrow energy intervals around 
E = Eq and E = E\, we consider resonant and non- 
resonant scattering separately. 

The paper is organized as follows. In Sec. II we derive 
a general expression for skew-scattering part of the 2D 
transport scattering cross section. This expression is lin- 
ear in spin-orbit coupling and contains all the scattering 
phases in the absence of spin-orbit coupling. In Sec. Ill 
we extract the energy dependence of the skew scattering 
in different domains for a model circular-barrier poten- 
tial. Resonant skew scattering is studied in detail in Sec. 
IV. In Sec. V. we consider the skew scattering in the 
Kondo regime. Concluding remarks are presented in Sec. 
VI. 



II. GENERAL FORMALISM 



As a result of the spin-orbit term 



H so — Act j^-L z 

P dp 



(2) 



in the Hamiltonian, where A is the spin-orbit constant 
and L z = —ihd/dO is the z-component of the orbital an- 
gular momentum, we have Si. a ^ <5;,-cr, where <5/ i(T are the 
scattering phases in the channel with orbital momentum 
I. 

Characteristics of the skew scattering, relevant for 
transport^ is 



I a (E) = 



:da a " 



(3) 



where da a /d6 is the asymmetric part of the differential 
scattering cross section 



da c {E) d<7 s (E) da a {E) 



d6 



dB 



dO 



(4) 



While ct is explicitly present in Eq. ([3]), I a (E) is, in fact, 
ct- independent. Expression for da a / 1 dO in terms of scat- 
tering phases reads 

= ^5> 24 ^ - l )^ 2lSi, ° i) - m (5) 

Eq. 



yields for I a (E) the following general relation: 
i sign(Z — I') 



(6) 



nk (cot Si yCr - i)(cot 5i>, a + i) 

= ^5Z{sin[2(5 !<tr - -, +1< J]-sm[2^]+sm[25, +1<t ,]}. 
i 

Obviously, the two last terms of the second line of Eq. ^ 
have vanishing contribution, and we will omit them later 
on. 

It is convenient to separate the spin-dependent and 
spin-independent parts of scattering phases 



5 La (E) = S°(E)+aSl(E). 



(7) 



Spin-orbit correction, Sl(E), can be expressed in terms 
of radial eigenfunctions of the continuous spectrum, 
Xi (p, E) , in the absence H so as follows 



51(E) 



2m\l 



d P ^ X Up,E). (8) 



P dp 

Large-p behavior of properly normalized Xl(Pi E) is 



X i(p,E) = cosSfWJwikp)- sin 5° (fyNwikp) 



(9) 
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Substitution of Eq. (O into Eq. ([6]) and expansion over 
Sf yields 



>2(«?-«?h). (10) 



Via simple transformation of the sum in Eq. (|10j) and 
insertion of the expression for 6} , Eq. ([5]) , we find 



I a {E) 



4toA 
irkh 2 



E 

/ — — oo 

,2 



dp dV(p) 
P dp 



(l + l) X f +1 (p,E) cosp^ -^)]. (11) 



In Eq. (fll|) , it is also convenient to make use of the iden- 
tity 



cos 2(5? - Sf^) = 1-2 



(tan<5° -tan^) 2 
(l + tan 2 < 5°)(l + tan 2 ( 50_ 1 )' 



(12) 



The advantage of using this identity is that the first term 
does not contribute to the sum in Eq. (TTTj) , and the ex- 
pression for I a (E) acquires the form 



dp dV(p) 



P dp 



lxf(p,E) 



(13) 



-(l-l)xli(p,E) 



(tan^-tan^) 2 
(1 + tan 2 <5°)(l + tan 2 5°~J' 



The form Eq. (fl3|) clearly illustrates that, for a weak 
scattering potential, skew scattering is oc V 3 . Indeed, 
when the scattering phases are small, the last fraction 
in Eq. (fT3")) can be replaced by (<5 ; — <5/_i) 2 oc V 2 , 
while for xi(p) one can use the free radial wave func- 
tions x1(P>E) = (np/2) 1/2 Ji(kp), Another power of V 



comes from 



dV(p) 



Eq. (fl3)) is also convenient for the 



analysis of the energy dependence of the skew scattering 
from a strong potential. We perform this study in the 
next Section. 



III. CIRCULAR-BARRIER POTENTIAL 

It is a textbook knowledge 20 that spin-independent 
scattering cross section from a strong potential, V ^S> 
h 2 jrab 2 , has different forms in three energy intervals: 
(i) low-energy scattering, E <C in which only zero 
orbital momentum, I = 0, contributes to the cross sec- 
tion; 

h 



(ii) intermediate energies (semiclassical regime), 

e V 



< 



vy (A) 



, where the scattering cross section is 
determined by high I 3> 1, and finally, 
(iii) high-energy scattering, E ^> V V/ {^gi^j 
corresponds to the Born approximation. 



which 



As we will see below, the skew scattering exhibits dif- 
ferent universal behaviors in the above three domains. 
We will see that, within the interval (ii), in contrast to 
the spin-independent scattering, the skew scattering has 
an additional scale at E ~ Vb, where it passes through 
the maximum value. Note that, for strong potential, 
the energy E ~ Vq is intermediate between h 2 /mb 2 and 
mV^b 2 1 'h 2 , which are the boundaries of the interval (ii). 

We will perform calculations for a model of circular- 
barrier potential V{p) = Vo9(b — p). In principle, 
this model allows to incorporate spin-orbit interaction 
nonperturbativelyi^. However, we will use the fact that 
this interaction is weak and evaluate Eq. (fTB")) expanded 
with respect to the coupling parameter, A. 

Within the model V(p) — V o 0(b — p), the phases 5f 
are found from matching at p = b the expressions for 
x'{p)/x{p) inside and outside the well. This yields 



JVi.(J*)/|,|(i*)-f/f, | (t*)JV|,|()*)- 



where 



v = yVo — k 2 , wo = 2mVo/H 2 . 



(15) 



It is convenient to rewrite Eq. (|14|) in a different form 
using the recurrence relations for Ji(z) and Ni(z) 



(A -S^) 



Also, within a model of a circular well, we have 4r- 

1 ' dp 

S(p — b), so that the general expression Eq. (JTTJ) assumes 
the form 



I a (E) 



4Xv 



(17) 



" (tan^-tan^) 2 
^ (1 + tan 2 <J°)(1 + tan 2 Sf_j) 



Here we have introduced a notation 
1 + tan 2 5f(k) 



I I -I 



Zi Zi-i 



[J w (kb)-t a n6f(k)N w (kb)Y 



(18) 



so that the spin-orbit corrections to the scattering phases 
are given by 



si 



ttXIvq 



2Z> 



(19) 



From Eq. HH), we can express Zu\ via the Bessel func- 
tions 



4 



Z\i\ - 



' ixkb 



kb k I\i\(vb) 



J\l\(kb) - J w+1 (kb) 



kb 



k Im(vb) 



N w (kb) ~ N w+1 (kb) 



(20) 



We start the analysis of Eq. ([17]) for the strong poten- 
tial Vo ^S> ^p- from the low-energy domain (i), kb <C 1. 
In this domain, spin-independent scattering is dominated 
by a single phase 



8° 



arctan 



21n(fc6) 



(21) 



which follows from z<Cl behavior of Nq(z). The phases 
with higher momenta are much smaller, namely 8\ ?(E) ~ 
(kb) 21 . On the other hand, parameters Zi grow rapidly 
with L 



v b- 



r 2 (0 



(22) 



as follows from Eq. (f20|) (here T(l) is the Gamma- 
function). For this reason we can replace the sum in 
Eq. ([T7]) by a single term, (<5°) jZ\, so that 

_ 4A« (<5 ) 2 _ 



I a (E) 



kb 
hi 2 (kb)' 



(23) 



First we emphasize that, unlike er s , the low-energy skew 
scattering is governed by two phases, 8® and 8^; the latter 
phase enters into Z\ , as can be seen from Eq. ([T8]) . This 
is a natural consequence of the form Eq. @ of the spin- 
orbit Hamiltonian. We also note that, compared to a s , 
which behaves as a S (E) oc l/ln 2 (fc6), Eq. ([23]) contains 
an additional factor (kb) oc E 1 / 2 , growing rapidly with E. 
Although derived for particular model of a circular-well 
potential, this low--E result is universal within a factor 
for a general "strong" potential, v$b 2 S> 1. 

1 /2 

Next we consider the energy domain 1 <C kb <C v b, 
which belongs to (ii). In this domain, the following sim- 



plifications are possible. Firstly, since vb 



1/2 



b » kb, 



we can replace I[(vb) / Ii(vb) in Eq. ([2"0"]) by 1. Such a re- 
placement is valid only for I <C vb. However, the relevant 
range of momenta in this domain is narrower, namely 
I < kb <C vb. This is because, for I < kb, the Bessel 
functions, Ji(kb) and Ni(kb) oscillate, while for I > kb 
they behave as Ji(kb) ~ (kb/l) 1 and N t (kb) ~ (kb/l)~ l . 
This behavior suggests that, tan^ falls off very rapidly 
once / exceeds kb, and the sum over I in Eq. (fTT)) should 
be terminated^! at / = kb. Note now, that in the interval 
1 < / < kb, the ratio l/(kb) in Eq. (fTB|) can be neglected 
compared to v/k. This yields the following simplified 
expression for the phases, Sf(E), 



-I, I , fv\ nl it 

tan oj = — tan { kb + arctan I — I 

\kJ 2 4 



(24) 



On the other hand, neglecting l/(kb) in Eq. ([20)) leads us 
to the conclusion that in the interval 1 < I < kb all Zi 
are equal to each other and are equal to 



Zi = 



Tkb 



irvob 
2k 



(25) 



Then the combination [l/Zi — (I — 1)/Z;_i] simplifies to 
2k/(TTVob). Lastly, we notice that 8f and S^_ x are related 



^-^ 1=7 r/2. 



(26) 



This follows from Eq. (|24|) . Using the latter relation, it 
is easy to see that all combinations containing phases in 
front of the square brackets in the sum Eq. ([TT)) are equal 
to 1. Thus, the sum over I in Eq. (|17[) is simply equal to 
2k 2 /(iTVo), and the result for I a (E) reads 



r(E) 



= ^kb. 

Kfc6<v^6 Trb 



(27) 



We see that the two asymptotes Eq. (J2"3"[) and Eq. (JS7J) 
match at kb ~ 1. In both limits I a (E) increases essen- 
tially as oc E 1 / 2 with E. Curiously, the magnitude of 
potential, vq, drops out from the asymmetric part of the 
scattering rate in both limits. 

The growth of I a (E) with energy is terminated at 
E ~ Vo, i.e., at kb ~ vb, which also belongs to do- 
main (ii). This can be seen from Eqs. (fT6|) . ([20]l in the 
following way. For I ~ kb ~ 3> 1 we can use the 

>l/2 



asymptotic expression |iJ(z)/J/(z)] « (l + l 2 /z 2 ^j 
the ratio I[(vb) / h(vb) . Then Eq. ((20]) yields 



irkb 
Zl = — 



Tb-^ 11 



2 b 2 



for 



(28) 



Using the same asymptotic form in the expression 
Eq. (fTH]) for phases, <5 Z ° amounts to replacement of v/k 
by (y/v 2 b 2 + I 2 — l) I (kb) in the argument of the arctan- 
gent in Eq. ([24]) . Still, the difference (8f — remains 
7r/2 with accuracy of a small parameter l/(kb). This 
is sufficient to replace by 1 all combinations containing 
phases in front of the square brackets in the sum Eq. (|17[) . 
Upon substituting Eq. ([28]) into Eq. (fTT]) , the summation 
over I can be easily performed, and we obtain 



I a (E) 



4Au 
k 
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FIG. 2: (Color online) Intermediate energy domain: 
h 2 /(mb 2 ) < E ~ Vo < mV£b 2 /h 2 . The skew-scattering 
contribution, I a (E), to the scattering rate, in the units 
SAy^/yr, is plotted from Eq. (|29[) versus dimensionless en- 
ergy x = E/V . 



width, r, of the quasilocal state. For the model poten- 
tial Fig. 1 the calculation of the width (with prefactor) 
is presented in the Appendix. Below we consider skew 
scattering near the resonance for two particular cases: 

1. Resonance for zero angular momentum, / = 0. In 
this case, the phase 5q(E) changes by tt as E is swept 
across the resonance. As we have seen above, in the 
low-energy limit, kb <C 1, the phases S® fall off rapidly 
with I. In calculating I a , the phases, S±i(E), should be 
retained, since there is no skew scattering without them. 
All phases 8i(E) with |2| > 2 can be neglected. 

2. Resonance for angular momenta, I = ±1. Now the 
phases S±i(E) exhibit resonant behavior. One has to 
retain the phase So(E); all phases Si(E) with |Z| > 2 can 
be neglected. 

Cases 1 and 2 are dramatically different because H so 
does not split the level / = 0, but does split levels I = ±1. 
For this reason, in the case 1, S] zl (E) can be calculated 
perturbatively from Eq. ©, so that for So we have 



So 



arctan 



r 



2(E - E ) 



(31) 



Now we see that I a {E) reaches a maximum at k = 

1 /2 

(vq/2) ' , which corresponds to E = Vq/2. The maxi- 
mal value is equal to 



I a (V /2) = 



4(V2 + 1) 



Ax/wo 



(30) 



In the case 2, So is non- resonant and is still not affected 
by H SOl while both scattering phases Si and #_x exhibit 
a resonant behavior 



S\, a = 5- 



1,-cr 



= — arctan 



2(E-Ex- SE h0 



(32) 



For smaller energies, Eq. (|2"9"|) reproduces the result 
Eq. flUI). Although Eq. ([29]) was derived for E < V , 
it remains applicable also for above-barrier scattering, 
E > Vq. It is seen that, after passing the maximum, 
I a (E) falls off as E' 1 / 2 . The energy dependence of I a {E) 
is illustrated in Fig. 2. 

We note that the position of maximum in I a (E) is 
model-dependent, in the sense, that for a general poten- 
tial with radius, ~ &, and magnitude, ~ Vo, the position 
of maximum can differ from Vq/2 by a numerical fac- 
tor ~ 1. However, the existence of maximum in I a (E), 
followed by oc E~ x l 2 decrease, is model-independent. 
Within a scaling factor, Eq. (f2"9"| applies up to the bound- 
ary, E ~ mVo 2 b 2 /ti 2 , of the domain (iii), where the Born 
approximation applies. From E ~ Vq to this boundary 
I a (E) drops by a large factor ~ \/vob 2 . In the high- 
energy tail (iii), the behavior of I a (E) depends strongly 
on the shape of the potential. 



Now we have to express the splitting SEi : ± a through H so . 



SEi <a = aSE\ = Xal 



J-dp^jRQ^E,)}' 
J^dppi^ME,)} 2 



(33) 



where i?i(/?, E) is the radial wave function of the discrete 
state, i.e., the coupling to continuum is neglected. Im- 
portant is that the width, of the state E\ is unaffected 
by the spin-orbit interaction. 

Now we evaluate Eq. ^ for the cases 1 and 2. Upon 
neglecting Si with \l\ > 2, we obtain from Eq.© a sim- 
plified form of I a (E) 

I a (E) = [sin(2(Si, ff ) + sin 2{S - Si, a ) (34) 
irk 

+ sin 2(5_i, ff -S )- sin(25_i )0 .)] 
4 

= — sinf^i — S-i) sin(i5o — Si — S-%) sin Jo- 
irk 



IV. I a (E) IN THE PRESENCE OF QUASILOCAL 
STATES 

The simplest model in which the quasilocal states 
emerge, is the repulsive scattering potential with attrac- 
tive core, as shown in Fig. 1. Quasilocal state does not 
affect the scattering process when the deviation of energy 
of the incident electron from the resonance exceeds the 



Since I a is independent of a, in the last identity we had 
dropped the subindex, cr, in the scattering phases. In 
the case 1, the phases S\ and S-\ are small, and Eq. 
assumes the form 



I a (E) = ^(S 1 - sin 2 So 

i r 2 

^^-^(E-Eoy+n/A- 



(35) 
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FIG. 3: (Color online) The shape of I" for small E, E 1 < 
h 2 /(mb 2 ), is plotted from Eq. (|3Tf) versus dimensionless devi- 
ation x = 2{E — E\)/Ti from the resonance, for three different 
values of dimensionless spin-orbit splitting, SEi/Fi. 



The difference (Si — <5_i) is proportional to spin-orbit 
constant, A, and can be evaluated from Eqs. (fT9| and 
(H2J), namely, 



Si -6- 



(36) 



Eq. (J35J) indicates that in the case 1 the energy depen- 
dence of the skew scattering is the same as the energy 
dependence of o~ s (E), i.e., it exhibits a resonant enhance- 
ment near E = E . 



In the case 2, So is small and nonresonant, while Si 
and (5_ i have resonances at E = E\ ± SEi . Upon using 
Eq. ([32]), we find 



I a (E) 



4Sq 
irk 



(37) 



T\(E-Ei)5Ei 



[(E -Ei- SEif + r?/4] [(E -Ei+ SEif + r?/4] " 



Our prime observation is that, in the case 2, the behaviors 
of a s (E) and I a (E) are vastly different. While a s (E) is 
oc T\/ [(E — Ei) 2 + T^/4] , the skew scattering is an odd 
function of (E—Ei). Actual shape of I a (E) is determined 
by the relation of two small energies SEi and Ti. For 
SEi <C Ti, position of a maximum (and minimum) and 
the width of I a (E) are all ~ T\. In the opposite limit, 
SEi ^> Ti, the extrema are located at (E — Ei) = ±SEi, 
while the width ~ Ti is much smaller. Behavior of the 
skew scattering near the resonance is illustrated in Fig. 3 
for both limits. 

At this point we emphasize a significant difference be- 
tween the two cases. In the case 1 we have I a tx dV/dp 
at p = b while in the case 2, the characteristic scale SEi 
is cx dV/dp at p = a. 

In the end of this Section we consider the case when 
the resonance E — Ei falls into the domain of "semi- 

1 /2 

classical" scattering, 1 <§; kb < v Q b. Nonresonant skew 
scattering is described by Eq. (|27|) in this domain and 
requires taking into account all the angular momenta up 
to I = kb. For the same reason, in calculating resonant 
skew scattering, one has to keep all four terms in the sum 
Eq. ([S]) which contain either Si or S-±. We have 



J 



61? (E) = \ [sin 2(6 1 -6 2 ) + sin 2(S - Si) + sin 2(<5_i - 6 ) + sin 2(<5_ 2 - <5_i 
nk 

2 

nk 



(38) 



= — [sin((5i - <5_i -6 2 + 6- 2 ) - s'm(Si - 6-i)] cos(<5i + 6-i - 2S ) = (5i 2 - Sl)[cos2(Si - S ) + cos2(<5_i - S )]. 

I 



In the second identity we used the fact that S 2 + 6- 2 ~ 
2Sq, which follows from Eq. (|24| : in the third identity 
we replaced sin(<52 — 6- 2 ), which is nonzero only due to 
spin-orbit-induced corrections, by (6 2 — <S^_ 2 ) ~ iXk/b. 
We note now, that the nonresonant parts of differences 
(Si — So) and (S-i — So) are 7r/2 and — 7r/2, respectively 
[see Eq. (|26p]. while the resonant parts of Si and <5_i are 
given by Eq. (|32|) . This leads us to the final result 



61? (E) « g 



(E - Ei- SEi) 2 -T\/A 



(E 



Ei-SEiY 
(E-Exl 



fTf/4 
SEif 



(39) 



rf/4 



(E-Ei + 5Eif + T\/A 



It is seen from Eq. (|39p that the characteristic mag- 
nitude of the resonant skew scattering contribution is 
8I?(E) ~ A/fe. On the other hand, the background value, 
I a (E) ~ Afc [see Eq. (|27)) ]. is much larger. This is because 
the background value is the sum of kb 3> 1 contribu- 
tions. Although 51? (E) constitutes a small correction, it 
has a lively energy dependence. This dependence is illus- 
trated in Fig. 4 for different values of dimensionless ratio 
5Ei/Ti. We see that, as 8E1/T1 increases, the structure 
in 81? (E) crosses over from one minimum to two minima. 
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FIG. 4: (Color online) The shape of the resonant contribution 
to the skew-scattering part of the scattering rate, 51? , in the 
domain, h 2 /(rnb 2 ) <£< Vo, is plotted from Eq. (|39|l versus 
dimensionless deviation, x — 2(E — Ei)/Ti, for three different 
values of dimensionless ratio SE\/Ti. 



V. MANY-BODY EFFECTS IN THE 
RESONANT SKEW SCATTERING 

As it was mentioned in the Introduction, the quantity 
relevant for transport is I a (E) at E = Ep. Then the 
resonance condition, Ep = E Q , can be satisfied cither 
if electron concentration or the potential of the attrac- 
tive core can be controlled. On the other hand, even 
away from resonance, at Ep > Eq, the symmetric part 
of scattering cross section, o- s (Ef), experiences a strong 
enhancement when the temperature, T, is lower than the 
Kondo temperature, Tk- This prominent many-body ef- 



fect stems from the Hubbard repulsion of two electrons 
in the quasilocal level. 

Spin-orbit interaction does not lift the degeneracy of 
the ground state, and thus, does not affect a s (Ep). On 
the other hand, since the skew scattering is essentially 
due to interference of I = and I — ±1 scattering chan- 
nels, one should expect the Kondo enhancement of a a 
alongside with enhancement of a s . 

Quantitatively, this enhancement can be found in the 
following way. Kondo effect modifies Eq. (|35|) . Namely, 
the scattering phase, So, should be replaced by the 
temperature-dependent Kondo phase, 8q(T). The low- 
T and high-T asymptotes of Sq(T) are the following 



1 - 



sin 2 So = 



(40) 



1 -• />/*. 



16 \n^(T/T K ) 



Eq. (|40|) describes how I a drops with increasing tem- 
perature from its "unitary" value I a = ir 2 \kp, where kp 
is the Fermi momentum, to the background value, given 
by Eq. flU) with k -> kp. 

A very lively behavior of the skew scattering emerges 
when the quasilocal level, E\, is close to the Fermi level, 
while Eq is well below the Fermi level. Then the Kondo 
resonance develops in I — channel at T < Tk, while the 
phases 5±\(Ep) are strongly sensitive to the deviation 
(Ei - Ep). 

A particularly interesting issue is what happens to the 
sign reversal of I a , found in the previous Section. In 
this case, both 5o and 5±i are resonant; sin<5o is defined 
by Kondo resonance Eq. (I40| . whereas S±i are given by 
Eq. (1211). Upon setting E — Ep in Eq. J32J) and substi- 
tuting it into Eq. (|M)) , we obtain 



I a (E u E F ,T) 



ASEtT! 
irk 



sin 2 <5 (T) 



Tx(Ep - E 1 )cot6 {T) + [(E F - E l - 8E{)(E F - E l +5E l ) - T 2 /4] 



[(E F -E x - SE-l) 2 
I 



T 2 /4] [(Ep-Ex+SExf + Tl/A] 



(41) 



At "high" temperature, 8q is small and we return to 
Eq. (|37|) . However, deep in the Kondo regime, the be- 
havior of I a (Ep) changes drastically. As it is seen from 
Eq. (j4Tj) . the dominant contribution to I a comes from the 
second term in the numerator, which is even function of 
(Ep — Ei), unlike the first term, which is odd. Therefore, 
upon decreasing temperature, the curve, I a (Ep), evolves 
from asymmetric to symmetric. 

For So = 7r/2, the shape of I a versus Ep, which is pro- 
portional to the electron density, is shown in Fig. 5. We 
see that symmetric shape undergoes a strong transfor- 
mation as the splitting SE\JT\ increases. For each value 



of splitting there are two points of the sign change posi- 
tioned symmetrically with respect to the point Ep = E\. 



VI. CONCLUDING REMARKS 

The Lorentzian shape of I a (E), described by Eq. (|35D . 
suggests that it should be accompanied by the Fano 
feature™, since I a (E) must assume its nonresonant value 
away from E = Eq. The issue of Fano resonance near 
E = Ei is more delicate. As can be seen from Eq. ((57)) 
at the Fano-resonance condition Si « £_i w Sq/2 the 
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Upon substitution 



0,2 




15 -10 


-b 
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10 


15 


sE/r, 

8E/P, 

SE/r, 


= 0.2 
= 1.5 
= 4 
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-0,8| 









FIG. 5: (Color online) The shape of I a in the Kondo regime, 
T <C Tk, and resonance with quasilocal level, E = Ex, is 
plotted from Eq. (|41fl versus dimensionless deviation a; = 
2(i5i — _Ef)/Fi of Ei from the Fermi level for three values 
of the dimensionless spin-orbit splitting, 5E\/T\, of the level 
E = Ej. 



Xi(p,E) 



Vp 1 

Eq. (|A1[) acquires a Hermitian form 

I 2 - 1/4 



X/ 



2m 



- U(p) 



(A2) 



Xi = 0. (A3) 



Consider now the potential, depicted in Fig. 1 with 
V(p) — for p > b. To calculate the width, T, we con- 
sider an auxiliary potential, V"(p), which coincides with 
V{p) for p < b and is a constant V(p) = Vq for p > b. 
Localized state, E = Ei in this potential is stationary. 
Denote with xi(p) the corresponding radial wave func- 
tion, so that 



2m 



Ei - V(p) 



I 2 



1/4 



Xi = 0. (A4) 



Upon multiplying Eq. (|A3|) by xi and Eq. (|A4|) by xu 
subtracting and integrating from a to b we obtain the 
relation 



factor sin(<5i — <5-i) turns to zero. This observation can 
be interpreted as "cancellation" of Fano resonances due 
to interference in I = 1 and I = — 1 channels. Simi- 
lar "cancellation" was pointed out in Refs. [HU In 
these papers the photo-current, caused by infrared ex- 
citation of electron either from impurity into the con- 
duction band 23 or between two Zeeman subbands^ of 
InSb in a strong magnetic field, was studied experimen- 
tally. In fact, there is a general similarity between the 
sine reversal of the skew scattering with energy and the 
sign reversal of photocurrent as a function of magnetic 
field, observed in Refs.|23,24, the underlying reason being 
the interference of a non-resonant and two split resonant 
channels. 



[xa'i 



XlXi 



2m (F 



Ei) dpxiXi- (A5) 



The fact that the difference (Ei — E{) is much smaller 
than Ei allows us to keep in the left-hand side only the 
contribution from p = b. Then it is convenient to rewrite 
Eq. (gSl as 



Xi Xi 



2m(Ei-Ei) 
h 2 Xi(b)xi(b) 



dpXi{p)Xi{p)- (A6) 



To find the imaginary part of Ei we use the continuity of 
the logarithmic derivatives of xi a t P = b 
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Xi 26 H+(kb)' 



(A7) 
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APPENDIX A 

In this Appendix we derive an analytical expression for 
the width of a 2D quasilocal state. The Schrodinger equa- 
tion for the radial part of the wave function, Ri(p, E), in 
a azimuthally-symmetric potential, V(p), reads 



where H^(kp) is the Hankel function corresponding to 
outgoing wave at p — > oo. Since xi{p) IS real, the second 
term in the left-hand side of Eq. (|A6|) does not contribute 
to ^smEi . The imaginary part of Ei originates from log- 
arithmic derivative of the Hankel function in Eq. (|A7|) . 
Final expression for T/ = Smfij emerges upon setting 
Xi(p) = Xi(p) = \/P K i(vp), where Ki is the Macdon- 
ald function, which is the solution of Eq. (|A4j) in the 
barrier region, and extending the upper limit of integra- 
tion in Eq. (|A6|) to infinity. Both steps are justified if 
^smEi <C Ei. Then we obtain 



„// 1 „/ ( 2m 

R " + -p^ + [-w 



Et-V(p) -L) Rl=0 , (Al) Ti(E) = 



K 2 {vb) 



irm Jf{kb)+N 2 {kb) ^dppK 2 {vp) 



(A8) 
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Further simplification is achieved when the core is wide 

1 11 

enough, so that va — [2m(Vb — Ei)a 2 ] /H is large. 
Then we can use the large-p asymptote of the Macdonald 
function in the integrand in Eq. (|A8|) , which yields 



In the opposite limit, va <C 1, the exponent in the width, 
T;, is exp(— 2vb), while the dependence on a enters into 
the prefactor. 



r«(£) 



2h 2 ve- 2v{ - h -^ 
~^rt> ' Jf{kb) + Nf{kb)' 



(A9) 
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